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We investigate possible tests of CPT invariance on the level of event rates at neutrino factories. 

We do not assume any specific model but phenomenological differences in the neutrino-antineutrino 
masses and mixing angles in a Lorentz invariance preserving context, such as it could be induced 
by physics beyond the Standard Model. We especially focus on the muon neutrino and antineutrino 
disappearance channels in order to obtain constraints on the neutrino-antineutrino mass and mixing 
angle differences; we found, for example, that the sensitivity \m 3 — m, 3 \ < 1.9 ■ 10~ 4 eV could be 
achieved. 

PACS numbers: 14.60.Pq 


I. INTRODUCTION 

The CPT theorem |lj is one of the milestones of lo¬ 
cal quantum field theory. It is based on such general 
principles as Lorentz invariance, the connection of spin 
and statistics, and the locality and hermiticity of the La- 
grangian. The SU(3) x SU(2) x U(l) Standard Model 
of Elementary Particle Physics (SM), for which the CPT 
theorem is valid, is in very good agreement with all ex¬ 
isting experimental data. Beyond the SM, like in string 
theory models or in models involving extra dimensions, 
CPT invariance could be violated [pj.pl. Thus, the search 
for possible effects of CPT violation is connected to the 
search for physics beyond the SM. Many different tests of 
CPT invariance have been carried out. So far, no CPT 
violation has been found and rather strong bounds on 
the corresponding parameters have been obtained Q. 

One of the basic consequences of the CPT theorem is 
the equality between the masses of particles and their 
corresponding antiparticles. A strong bound on a possi¬ 
ble violation of CPT invariance has been obtained from 
the K°-K° system. This violation is characterized by the 
parameter 

_ TL R o. R o - ,... 

~ 2(A l - A S ) ’ U 

which can be related to measurable quantities ||. In 
Eq. (@), A L ,s = mL,s ~ ^l,s, m L ,S and T L>S are 
the masses and the total decay widths of the and 
K°s mesons, respectively, and Ti is the effective non- 
Hermitian Hamiltonian of the K°-K° system in the 
representation |AT°) and |A'°), which are eigenstates of 
the Hamiltonian of strong and electromagnetic interac¬ 
tions. For the complex diagonal matrix elements, we 
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have Ti R a. R o — 0 i^r^o and TC R o. R o — m R o ^T/yo, 
where m K o R o and T K o R o are the bare masses and the 
total decay widths of the K° and K° mesons, respec¬ 
tively, with corrections due to weak interactions. The 
CPLEAR experiment obtained [Q 

\m K o -m R0 | = (-1.5 ±2) • l(T 18 GeV. 


Using all relevant data on the K°-K° system, it follows 
that Q] 

\m K o - m R o| < 10 _i8 

^average 

where m average = {fn R o + fo)/2. Recently, also an 

upper bound on the mass difference between the B J] and 
B° mesons has been obtained [jTj] 


lm B o - m B o\ 
m B o 


< 1 . 6 - 10 " 14 . 


Here, we will consider possible CPT invariance tests 
that can be performed in future high-precision experi¬ 
ments with neutrinos from neutrino factories, which are 
now under active investigation j|, 0, [T7j. |TT|. |u| [b5| |T^, 
|Tt| . They have mainly been proposed to study neutrino 
oscillations in detail. In addition, in the framework of 
Lorentz non-invariant models, possible CPT invariance 
tests with neutrino experiments have been discussed (see, 
e.g., Refs. fu| [H]]). 

Compelling evidence for neutrino oscillations has been 
found by atmospheric Q and solar |lj 20, 0 neu¬ 
trino experiments. The following best-fit value for the 
atmospheric mass squared difference Am atm has been ob¬ 
tained 1513 ]: 


Ami ^ 2.5 ■ 10 


-3 


eV 2 . 


From the global analysis of all solar neutrino data, sev¬ 
eral allowed regions in the neutrino oscillation parameter 
space have been found. For the preferred so-called large 
mixing angle (LMA) solution in Ref. |^2| , the solar mass 
squared difference has been determined to be 


A m% ~ 4.5 ■ 10“ 5 eV 2 . 
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Furthermore, there are at present indications for neutrino 
oscillations with an even larger mass squared difference, 
which were found by the LSND experiment |2^]. From 
the analysis of the data of the LSND experiment, the 
best-fit value of the neutrino mass squared difference [Q 

A to lsnd — 0.24 eV 2 


was found. 

The strongest kinematical bound on the absolute neu¬ 
trino mass scale m\ is obtained from the endpoint of 
the /3-spectrum of 3 FI. The latest measurements yielded 
mi < 2.2 eV |25|, |b). From neutrinoless double /3-decay 
there exists also a strong bound |(m)| = \ — 

(0.2 — 0.6) eV for Majorana masses (for an overview see, 
e.g., Ref. ^7j). Here U e i are matrix elements of the neu¬ 
trino mixing matrix U and m; are the masses of the neu¬ 
trino mass eigenstates. Furthermore, somewhat weaker 
but similar bounds emerge from astrophysics and cosmol¬ 
ogy. It nevertheless follows from the existing neutrino 
data that neutrino masses are not equal to zero and that 
they are much smaller than the masses of all other fun¬ 
damental fermions (leptons and quarks). From empirical 
lepton and quark mass patterns a hierarchical (or inverse 
hierarchical) mass pattern seems to be rather plausible 


It is a general belief that the smallness of the neutrino 
masses requires some new mechanism beyond the SM. 
The classical mechanism of neutrino mass generation is 
the see-saw mechanism [^9|, which connects the small¬ 
ness of the neutrino masses with the violation of lepton 
numbers at an energy scale much higher than the elec- 
troweak scale. In this case, massive neutrinos have to be 
Majorana particles and the neutrino masses have to sat¬ 
isfy a hierarchy relation. The see-saw mechanism is based 
on local quantum field theory, and therefore, violation of 
CPT invariance cannot be expected. 

Furthermore, it has recently been suggested |3(J that 
the smallness of the neutrino masses could have a natu¬ 
ral explanation in models with large extra spatial dimen¬ 
sions. In such models, the smallness of the Dirac neu¬ 
trino masses follows from the suppression of Yukawa in¬ 
teractions of the left-handed neutrino fields, localized on 
a three-dimensional brane, and the singlet right-handed 
neutrino fields propagating together with the gravita¬ 
tional field in a bulk. In models with n extra dimensions, 
the neutrino masses are proportional to 


ri _ m „ 1Q ~i6 m 
V M n V n M g TeV’ 

where V n is the volume of the extra space, Mq — 
1.2 • 10 19 GeV is the Planck mass, and M ~ 1 TeV is the 
Planck mass in the 4 + n dimensional space. Moreover, 
there are other approaches to the generation of small 
Dirac or Majorana neutrino masses in models with extra 
dimensions (see, e.g., Refs. |n], ^|). Since the symme¬ 
tries of the SM are violated in the bulk, neutrino mass 
generation in extra dimension models is a plausible can¬ 
didate for the violation of CPT invariance p3[. 


In order to accommodate all existing neutrino oscilla¬ 
tion data, including the data of the LSND experiment, 
it is necessary to have three independent mass squared 
differences. Thus, we need to assume that there exist (at 
least) four massive mixed neutrinos, i.e., in addition to 
the three active flavors is e , u^, and i/ T at least one sterile 
neutrino has to exist Q. 

In Refs. H|, it was assumed that CPT violation in 
the neutrino sector can be so strong that the mass spec¬ 
tra of neutrinos Vi and antineutrinos T\ are completely 
different. In this case, it is possible to describe atmo¬ 
spheric, solar, and LSND neutrino data with a framework 
of three massive neutrinos and three massive antineutri¬ 
nos (assuming that AmL SND belongs to the antineutrino 
spectrum). Such an extreme picture can, in principle, be 
tested by the future MiniBooNE |h|, KamLAND [p7| , 
and other similar neutrino experiments [f33| . 

In Ref. JItJ , the effect of a term in the neutrino Hamil¬ 
tonian violating CPT and Lorentz invariance has been 
considered and the v lt —> and Pjf —* Pjf transition 

probabilities with the and T7^ coming from neutrino 
factories have been calculated. It was demonstrated that 
in such a model the effects of CPT violation could be 
rather large in a wide range of the corresponding param¬ 
eter values. 


II. BASIC FORMALISM 

In this paper, we will assume Lorentz invariance and 
consider possible violation of CPT invariance by the 
mechanism of neutrino mass generation. In the case of 
the usual neutrino mixing, we have 

^olL ^ v U a i l / iL 1 (2) 

i 

where U is a unitary mixing matrix and ig are the neu¬ 
trino fields (Dirac or Majorana) with masses rrii. The 
neutrino flavor state \v a ) is given by 

K) = \vi\muL), (3) 


where |z/j; nii , L) are the neutrino states with masses TOj, 
negative helicity L, 3-momentum p, and energy 

Ei = \J mf + p 2 ~ p + ^ (4) 

in the ultra-relativistic limit. |f8[ For the antineutrino 
flavor state \V 1 a) we have 

\T^) = ^2u ai \W,m i: R) (5) 

i 

in the case of Dirac neutrinos and 

[Ic) = \vj-,mi,R), (6) 
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in the case of Majorana neutrinos. In these relations, 
\Ui]nrii,R) and \iSi;mi,R) are the right-handed antineu¬ 
trino and Majorana neutrino states, respectively, which 
also have the 3-momentum p and the energy Ei. 

Assuming the usual Lorentz invariant propagation of 
neutrino states for the neutrino and antineutrino transi¬ 
tion probabilities in vacuum, we find the expressions 


P(r 


') = 


E^ 


JJ* 


( 7 ) 


and 


') = 


E 


77*,. e 

ol'i ° 


-iA mf, 7TE7 


tl 2JS TT„ 
v ry 


( 8 ) 


which automatically satisfy the relation 


P(l/ a > L , a f ') — P(t / a / * Z-'a). (9) 

In Eqs. (f7|) and (||), Am?- = m 2 — m 2 is the mass squared 
difference, L ~ t is the distance between the source and 
detector, and E is the neutrino energy. Note that Eq. (|[) 
is a consequence of CPT invariance inherent to standard 
neutrino mixing and oscillations. 

If the generation mechanism of neutrino masses and 
mixings violates CPT invariance, then the relations for 
antineutrino flavor states will differ from Eqs. (g) and (g). 
In the case of massive Dirac neutrinos, the antineutrino 
masses rrii will be different from the neutrino masses rrii , 
and the mixing matrices will, in general, not be connected 
by complex conjugation. Thus, for the antineutrino fla¬ 
vor states we have 


| i^ot) — ^ ^ U a i l^i, rrii, R). (19) 

i 


In the case of massive Majorana neutrinos, neutrinos 
and antineutrinos are identical. For the right-handed an¬ 
tineutrino flavor states, we therefore have 


7 )=E^ 


ii , R). 


( 11 ) 


Further on, we will assume that there is no violation of 
Lorentz invariance in the propagation of massive neutri¬ 
nos and antineutrinos. 


III. CPT TESTS AT NEUTRINO FACTORIES 

In this section, we will investigate the sensitivity of 
future high-precision neutrino oscillation experiments 
at neutrino factories to neutrino-antineutrino mass and 
mixing angle differences. Neutrino factories [g g will 
allow to investigate the phenomenon of neutrino oscilla¬ 
tions, which has been observed by the atmospheric and 
solar neutrino experiments, with unprecedented accu¬ 
racy. It will be possible to determine the leading neutrino 


oscillation parameters A?n| 2 and sin 2 2023 governing the 
z/p —> v T oscillations in the atmospheric region very well. 
Depending on their values, it will also be possible to limit 
or to measure the mixing angle $13 to search for the con¬ 
nected matter effects and to discriminate between a hi¬ 
erarchical neutrino mass spectrum and a mass spectrum 
with reversed hierarchy. In the most likely LMA case, 
the effects of CP violation in the lepton sector can be 
studied. Details of neutrino factory phenomenology can 
be found in Refs. JTo], 11, [TJ. [Tg [Tj. [Tg. As we will 
show below, because of the high precision of neutrino 
factories, we can estimate the sensitivity of experiments 
to the presumably small violations of CPT invariance in 
the neutrino sector, being an unambiguous sign of new 
physics. 

At neutrino factories neutrinos will be produced in 
muon decays fi + —> e + v e T7^ (or fi~ —> e~T7^v ll ). The 
straightforward way to test CPT invariance at neu¬ 
trino factories would be to check the appearance relation 


P(r 


J = P(7 


(or P(i 


:) = p(* 


u e )) with neutrinos from fj, + (/r _ ) decays. However, such 
tests would require to measure the sign of the charge 
of the produced lepton. The sign of a muon charge can 
be determined very reliably, but measuring the sign of an 
electron (or positron) charge is a rather challenging prob¬ 
lem. The possibility to measure the electron (or positron) 
charge with moderate efficiency with liquid argon detec¬ 
tors would not be precise enough. Therefore, we consider 
a CPT invariance test in the v fl and TjJ disappearance 
channels by checking the equality 


P K -> *v) = p (^ -*• *V)- 

The and Pjf disappearance channels have several ad¬ 
vantages: 

1. The effect of neutrino oscillations in the atmo¬ 
spheric mass squared difference region is large. 

2. The matter effects are small. 


3. There is no relevant background from the T^T’s 
(v e ’s), which are accompanying the s (FjTs) in 
the decays of the /u - ’s (/x + ’s). 

4. The event rates are high for obtaining good statis¬ 
tical information. 

We will only consider the possible violation of CPT in¬ 
variance in the v tJ —> and Pjf —*► TjJ oscillations. If CPT 
invariance is violated, then these oscillations will be char¬ 
acterized by the leading parameters Am ^ 2 , sin 2 2023 and 
Am| 2 , sin 2 20 2 3 , respectively. 

In Ref. [HH, a comprehensive study of the accuracy 
of the measurement of neutrino oscillation parameters in 
neutrino factory experiments was performed. Our calcu¬ 
lations will be based on this study. Since matter effects 
give only small contributions to the v IL and Pjf survival 
probabilities, uncertainties in the Earth matter density 
profile are of little importance for the parameter mea¬ 
surements. In Ref. |H|, Fig. 3, the relative statistical 
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errors of the parameters 5 Ato 32 and <56*23 7 determined by 
a general analysis including correlations, are plotted as 
functions of the luminosity 

C = 2 Nfj 777-kt, 

where N /t is the number of stored muons per year and 
Wkt is the mass of the detector in kilotons. 

Violation of CPT invariance in neutrino oscillations 
can be characterized by the following parameters: 

<5 = | A 77132 — Ato 32 |, (12) 

e = | sin 2 2023 - sin 2 20231 ■ (13) 

If the minimal neutrino mass mi and the CPT violat¬ 
ing effects are small (mi <C Ajti.q , jm 3 — m 3 | <C 

(7U 3 )average); then we find for the hierarchical neutrino 
mass spectrum or the spectrum with reversed hierarchy 
that 


<5 ~ 2 acpT A?n 32 , 


where 


acpT 


m 3 - m 3 | 


(m 3 ) 


average 


(14) 


(15) 


is a dimensionless parameter which characterizes the vi¬ 
olation of CPT invariance. We can also write e as 


e ~ 


2 &cpt \J sin 2 20 23 ^/l — sin 2 20 : 
= 2 &CPT #23 Sin 4023, 


'23 arcsm \J sin 20 2 3 
(16) 


where 


#CPT = 


I #23 — #23! 

(# 23 ) average 


(17) 


The experimental sensitivity to the possible CPT vio¬ 
lation is given by the accuracy with which the param¬ 
eters acpT and/or &cpt can be measured. In order to 
estimate the sensitivity we will treat the neutrino and 
antineutrino channels as different experiments which are 
not combined to fit a common Am§ 2 and 023- In order to 
establish an effect we therefore need to compare the val¬ 
ues of the parameters acpT and 6cpt, which are describ¬ 
ing the asymmetry between these two experiments, with 
the corresponding statistical errors of the neutrino oscil¬ 
lation parameters determined in Ref. [p~5| , Fig. 3. Only 
if the mass squared or mixing angle difference between 
neutrinos and antineutrinos is larger than the respective 
relative statistical error 5 Ato§ 2 or 6623 of the measure¬ 
ment of Am 32 or 023, CPT violation will be detectable 
on the respective confidence level of the statistical evalu¬ 
ation, i.e., the sensitivities <5ocpt and Sbcpr to the CPT 
violating parameters ocpt and &cpt are given by: 


<5acpT 

Sbc.PT 


5Ato| 2 
2 ! 
5023, 


(18a) 

(18b) 


where ocpt < 5ocpt and &cpt < 56cpt- The factor of 
two in the first relation comes from the translation from 
mass squared differences to masses for a hierarchical (or 
inverse hierarchical) mass spectrum in Eq. ©• As an 
example, a statistical error of 7% in the determination of 
Am 32 would correspond to a mass asymmetry sensitivity 
between neutrinos and anti neutr inos o f 3.5 %. The sen¬ 
sitivities described by Eqs. (18a) and (18b) are plotted 
in Fig. [I], where the sensitivity 5 ocpt to the asymmetry 


FIG. 1: The sensitivitities SacpT and <5&cpt of an estimate 
of the asymmetries ocpt and &cpt at a neutrino factory as 
functions of the luminosity £. The solid curve refers to the 
mass asymmetry ocpt (hierarchical or inverse hierarchical 
mass spectrum only) and the dashed curve to the mixing an¬ 
gle asymmetry &cpt • The underlying calculations in Ref. |jl5[| , 
Fig. 3, were performed with 50 GeV muon energy and base¬ 
lines of 7000km (023) and 3000 km (Am§ 2 ). 



ocpt is shown as a function of the luminosity C by the 
solid curve. The respective statistical errors were calcu¬ 
lated for a muon energy of E = 50 GeV and for a base¬ 
line length of L = 3000 km. Similarly, the dashed curve 
shows the sensitivity 5&cpt to the asymmetry &cpt, with 
the statistical errors calculated for E = 50 GeV and 
L = 7000 km. From these curves, one can, for example, 
read off for a 10 kt detector and 10 20 stored muons per 
year during 5 years (c/., the vertical line in the plot) that 
ocpt ^ 3.8 • 10 -3 and &cpt 4.3 • 10 -2 . For the mass 
difference of neutrino and antineutrino we then obtain for 
a hierarchical or inverse hierarchical neutrino mass spec¬ 
trum 1 777.3 - m 3 1 = acPT(TO3)average — acPTA/Am 2 tm < 

1.9 • 10” 4 eV. 


IV. SUMMARY AND CONCLUSIONS 

CPT is a fundamental symmetry preserved in any 
Lorentz invariant local quantum field theory. Especially, 
the SM is a CPT invariant theory. However, CPT in- 
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variance can be violated in models beyond the SM, like 
models with extra dimensions or string theory models. 
It is important to note that the expected effects of CPT 
violation depend on the assumed model. If the Planck 
mass is close to the TeV scale, such as it is for models 
with large extra dimensions, these effects could be ob¬ 
servable in future experiments. We especially addressed 
the question of CPT violation by small neutrino mass or 
mixing angle differences between neutrinos and antineu¬ 
trinos, which could, most plausibly, be generated by a 
mechanism beyond the SM. Furthermore, we investigated 
the sensitivity of future neutrino factory experiments to 
the presumably small mass and mixing angle differences 
in the and TTJf disappearance channels. Finally, we 
have shown that for the neutrino-antineutrino mass dif¬ 
ference in a hierarchical (or inverse hierarchical) neutrino 


mass spectrum, the upper bound 

|m 3 — m 3 | < 1.9 • 10 -4 eV 
can be obtained. 
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